arXivil504.04134vl [cond-mat.mtrl-sci] 16 Apr 2015 


Dipole-Allowed Direct Band Gap Silicon Superlattices 


Young Jun Oh, Sunghyun Kim, and K. J. Chan^ 

Department of Physies, Korea Advaneed Institute of Seienee and Teehnology, Daejeon 305-338, Korea 

In-Ho Lee 

Korea Researeh Institute of Standards and Seienee, Daejeon 305-340, Korea and 
Center for In Silieo Protein Seienee, Sehool of Computational Seienee, 

Korea Institute for Advaneed Study, Seoul 130-722, Korea 

Jooyoung 

Center for In Silieo Protein Seienee, Sehool of Computational Seienee, 

Korea Institute for Advaneed Study, Seoul 130-722, Korea 
(Dated: April 17, 2015) 

Silicon is the most popular material used in electronic devices. However, its poor optical properties 
owing to its indirect band gap nature limit its usage in optoelectronic devices. Here we present the 
discovery of super-stable pure-silicon superlattice structures that can serve as promising materials 
for solar cell applications and can lead to the realization of pure Si-based optoelectronic devices. The 
structures are almost identical to that of bulk Si except that defective layers are intercalated in the 
diamond lattice. The super lattices exhibit dipole-allowed direct band gaps as well as indirect band 
gaps, providing ideal conditions for the investigation of a direct-to-indirect band gap transition. 
The transition can be understood in terms of a novel conduction band originating from defective 
layers, an overlap between the valence- and conduction-band edge states at the interface layers, and 
zone folding with quantum conhnement effects on the conduction band of non-defective bulk-like 
Si. The fact that almost all structural portions of the superlattices originate from bulk Si warrants 
their stability and good lattice matching with bulk Si. Through first-principles molecular dynamics 
simulations, we confirmed their thermal stability and propose a possible method to synthesize the 
defective layer through wafer bonding. 


I. INTRODUCTION 

Silicon is an important element used in modern elec¬ 
tronic devices owing to its abundance, feasibility for 
large-scale fabrication, easy formation of native oxide, 
and doping controllability of both electrons and holes. 
However, the optical property of Si is rather poor owing 
to its indirect electronic band gap nature. In an indirect- 
band-gap material like cubic-diamond Si (denoted as c- 
Si), optical transitions at the threshold energy occur only 
via momentum-conserving phonons. Therefore, the solar 
spectrum pertaining to the energy below the direct band 
gap of c-Si, z.e., approximately 3.4 eV, cannot be effec¬ 
tively absorbed without phonon assistance. 

To improve the optical property of Si, considerable 
efforts have been made; for instance, by introducing 
defects, such as erbium atoms, dislocations, and grain 
boundaries as recombination centers m, and/or by en¬ 
gineering the electronic band structure through nanopat- 
terning HE], nanostructuring [6], alloying with group- 
IV elements H Ej, applying strain, or combinations of 
these HEg. Very recently, several Si crystals with di¬ 
rect and quasidirect band gaps were computationally de¬ 
signed [T 3 HT 5 ] and their optical absorption properties 
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were shown to be significantly improved compared with 
c-Si: quasidirect band gaps are defined as Eg < Eg < 
Eg P 0.1b eY following a previous study [13], where Eg 
and Eg denote direct and indirect band gaps, respec¬ 
tively. However, these metastable structures were of rel¬ 
atively high energies, ranging 0.1-0.3 eV per atom, be¬ 
cause of distorted tetrahedral bonds, and the experimen¬ 
tal realization of their synthesis has yet to be explored. 


In this work, using a computational search method, we 
find pure Si-based superlattices that exhibit direct and 
optically allowed band gaps with good lattice matching 
with c-Si. The superlattice structure is composed of al¬ 
ternating stacks of bulk-like Si(lll) layers and a defec¬ 
tive layer containing Seiwatz chains [16]. The electronic 
structure evolves as the bulk-like Si portion increases, 
exhibiting a transition from the direct to indirect band 
gap. In superlattices with direct band gaps, the optical 
transition at the threshold energy is greatly enhanced. 
Hence we suggest that these superlattices could be used 
in thin-film solar cells. In the next section, we explain the 
details of the computational method. Then we present 
the structural, electronic, and optical properties of our 
designed superlattices. Finally, we discuss the thermal 
stability and the possible synthesis through wafer bond¬ 
ing, based on molecular dynamics simulations. 
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FIG. 1. Superlattices composed of Si(lll) layers (green) and a defective layer of Seiwatz chains (orange) are shown for (a) n 
= 3, (b) arbitrary n, and (c) n = 1. The black parallelepiped represents the unit cell spanned by the lattice vectors, ai, a 2 , 
and as. (d) The orientation of Seiwatz chains on the Si(lll) 2x1 surface is shown. 


II. CALCULATION METHOD 

We explored Si crystal structures with optically ac¬ 
tive direct band gaps by using the computational search 
method in which electronic properties are initially 
assigned and target materials are subsequently searched. 
We employed a combined approach m of conforma¬ 
tional space annealing (CSA) [THHH] for global opti¬ 
mization and first-principles calculations within the den¬ 
sity functional theory (DFT) framework. Without us¬ 
ing any specific knowledge of known Si crystals, we op¬ 
timized the degrees of freedom including atomic posi¬ 
tions |^/| and six lattice parameters (a, 6 , c, a, and 

7 ). The objective function used for selection in CSA is 
—+max [ 0 , {Eg — F^p], which promotes the formation 
of direct band gap. For each conformation, the enthalpy 
was minimized by performing DFT calculations which 
used the functional form of Perdew, Burke, and Ernz- 
erhof (PBE) [22] for the exchange-correlation potential 
and the projector augmented wave potentials [23], as im¬ 
plemented in the VASP code [24]. The wave functions 
were expanded in plane waves with an energy cutoff of 
400 eV. With a /c-point mesh using the grid spacing of 
X 0.02 A ^, the crystal structures were minimized until 
all forces and stress tensors were less than 0.01 eV/A and 


1.5 kbar, respectively. Einally, a twice finer A:-point mesh 
was used to determine the nature of band gap (direct vs 
indirect). 

The band gap sizes of semiconductors and insulators 
are usually underestimated with the PBE functional. Eor 
more accurate calculations for the band gaps and absorp¬ 
tion coefficients, we additionally performed quasiparti¬ 
cle calculations in the Go Wq approximation [25l |26] and 
solved the Bethe-Salpeter equation considering up to 
12 occupied and 16 unoccupied bands around the Eermi 
level, which were shown to be sufficient to ensure the 
numerical convergence Ha¬ 
ul. RESULTS AND DISCUSSION 

After an extensive search for Si crystal structures with 
optically active direct band gaps, we obtained a very dis¬ 
tinctive superlattice structure, especially for the system 
containing 14 Si atoms per unit cell. The structure con¬ 
sists of alternating stacks of three Si(lll) layers and a 
defective layer along the [ 111 ] direction of c-Si [Eig.j^a)]. 
The defective layer contains the so-called Seiwatz chains 
(SCs) [16], which are known to be formed via 2 x 1 recon¬ 
struction on the Si(lll) surface with half-monolayer cov¬ 
erage. The SCs lead to open channels consisting of five- 
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TABLE L For the Si(lll)n/Si(SC) superlattices with the cubic- and hexagonal-stacking sequences of the Si(lll) layers, the 
lattice type, the number of atoms per unit cell (iV), the energy relative to cubic-diamond Si (E), the type of band gap, the direct 
band gap size (Eg), and the indirect band gap size (EJ) are shown, based on the PBE calculations. The quasiparticle GqWq 
gaps are also shown for comparison. Here D and QD in parentheses denote direct and quasidirect band gaps, respectively, and 
lattice types are abbreviated, such as SM: simple monoclinic, SO: simple orthorhombic, and BCO: base-centered orthorhombic. 


n 

lattice 

N 

E (meV/atom) 

PBE-S(| (eV) 

PBE-Ej (eV) 

Go Wo-Ei, (eV) 

Go Wo-Ei (eV) 

1 

BCO 

6 

89 (QD) 

0.431 

0.430 

0.894 

0.847 





Cubic-diamond stacking 



2 

SM 

10 

46 (QD) 

0.906 

0.869 

1.346 

1.316 

3 

SM 

14 

42 (D) 

0.807 


1.197 


4 

SM 

18 

32 (D) 

0.832 


1.283 


5 

SM 

22 

26 (D) 

0.782 


1.218 


6 

SM 

26 

22 (QD) 

0.788 

0.774 

1.224 

1.215 

7 

SM 

30 

19 (QD) 

0.761 

0.741 

1.198 

1.185 

8 

SM 

34 

17 (QD) 

0.746 

0.726 

1.185 

1.166 

9 

SM 

38 

16 (QD) 

0.738 

0.712 

1.177 

1.154 

10 

SM 

42 

13 (QD) 

0.725 

0.699 

1.167 

1.147 





Hexagonal-diamond stacking 



2 

SO 

10 

72 (QD) 

0.562 

0.527 

0.980 

0.961 

3 

SM 

14 

49 (D) 

0.615 


1.049 


4 

SO 

18 

30 (D) 

0.578 


1.008 


5 

SM 

22 

34 (D) 

0.497 


0.914 



and eight-membered rings in the defective region. Note 
that all Si atoms are four-fold coordinated, without any 
coordination defects, and nearly ideal tetrahedral bonds 
are formed in the non-defective region. Owing to the 
almost identical structural match between the Si(lll) 
layers and c-Si, many additional superlattice structures 
[denoted as Si(lll)n/Si(SC)] can be constructed by vary¬ 
ing n, as shown in Fig. Bb)- For n = 1, the crystal 
is entirely composed of five- and eight-membered rings 
[Fig. B c)], very similar to an orthorhombic allotrope of 
Si, Si 24 , which has been recently synthesized [28] . 

In the cubic-diamond (Si-I) phase, the stacking se¬ 
quence of the Si(lll) layers is ABC ABC .while that 
of the hexagonal-diamond (lonsdaleite, Si-IV) phase is 
ABABAB.... In our superlattice system, various stack¬ 
ing sequences including the cubic and hexagonal ones are 
possible, as in SiC polytypes. In addition, in each de¬ 
fective layer, two configurations of SCs exist, which are 
related to each other by the translation of (ai + a 2 )/ 2 , 
where ai and ai are the lattice vectors on the 2x1 basal 
plane [Fig. [^d)]. Thus, depending on the relative po¬ 
sitions of two adjacent defective layers, two Bravais lat¬ 
tices can be formed, simple monoclinic (SM) and base- 
centered monoclinic (BCM) for n > 2, with some ex¬ 
ceptions for the case of hexagonal stacking (Table . By 
examining various superlattices, we found that BCM su¬ 
perlattices tend to yield quasidirect/indirect band gaps, 
whereas the nature of the band gaps (direct vs quasidi¬ 
rect) in the SM ones varies with n. 

Because c-Si is the most stable form of Si, we fo¬ 
cused on SM superlattices with cubic-diamond stacking; 
throughout this paper, unless otherwise specified, the cu¬ 
bic stacking is considered. For superlattices with n up 


to n = 13, the valence band maximum (VBM) is al¬ 
ways located at the F point, the center of the Brillouin 
zone (BZ). For n = 3-5, we found direct band gaps at 
F (Figs. i and [^, which were estimated to be 0.807, 
0.832, and 0.782 eV, respectively, using the PBE func¬ 
tional for the exchange-correlation potential (Table . 
As n increased, a direct-to-quasidirect band gap transi¬ 
tion occurred. The energy differences between the di¬ 
rect and indirect band gaps were smaller than 30 meV 
for 6 < n < 10. Although the band gap size decreased 
with increasing n, it was larger than that (0.62 eV) of 
c-Si, owing to the quantum confinement effect. When 
more rigorous quasiparticle Go Wq calculations were per¬ 
formed, the nature of the band gap did not change. The 
Go Wo band gaps were 1.15-1.28 eV (n = 3-10), close to 
the optimal value (1.1 or 1.3 eV) for solar cell applica¬ 
tions [29] (Table [^. For comparison, we also examined 
SM superlattices with hexagonal-diamond stacking and 
found direct band gaps at F for n = 3-5. However, the 
band gap sizes were reduced by approximately 0.19-0.29 
eV (Table |I|, similar to SiC polytypes, where the band 
gap tends to decrease as the number of hexagonal layers 
in stacking sequence increases. 

In our superlattice system, the lowest conduction band 
exhibits nearly flat dispersion along the F-X and F-Y 
directions (Figs.j^an d§. To understand the origin of the 
fiat dispersion around F, we considered the defective and 
Si(lll) layers separately (Fig.|^. In a slab geometry with 
a single defective layer sandwiched between two Si(lll) 
layers, a vacuum region was included in the same unit 
cell and Si dangling bonds were passivated by hydrogen. 
It is clear that the fiat dispersion is attributed to the SCs 
in the defective layer, as shown in Fig. [^ On the other 
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FIG. 2. (a) Symmetry points and lines in the Brillouin 

zone of the simple monoclinic lattice are shown, (b) The 
PBE band structure of the Si(lll)n=5/Si(SC) superlattice is 
shown. The thickness of red colored bands represents the 
degree of confinement in the defective layer, indicating that 
the lowest conduction band along the X-T-Y-C line is mainly 
derived from the Seiwatz chains. For the highest valence and 
lowest conduction bands along the X-F-Y-C line, the size of 
blue circles is proportional to the degree of contribution from 
the interface layer. 

hand, the energy states derived from the Si(lll) layers 
are affected by the zone folding and quantum confinement 
effects. In c-Si, the conduction band minimum (CBM) is 
located at six A-valleys close to the X-points in the BZ 
of the face-centered cubic lattice. When a SM cell with 
a 2x1 basal plane is adopted, two of the six X-points are 
folded to F in the monoclinic BZ. In contrast to c-Si, the 
hexagonal symmetry in the superlattice is perturbed by 
the SCs [Fig. [^d)]; hence, six A-valley states split into 
two folded states, Ai and A 2 , 3 , which have two- and four¬ 
fold degeneracies, respectively. The folded Ai points are 
located near F, whereas the folded A 2,3 points are close 
to A or E in BZ, depending on n [Fig. [^a)]. The lowest 
conduction bands at F, Ai, and A 2,3 are shown as a 
function of n in Fig. ^ While the A 2,3 states are mainly 
confined in the Si(lll) layers, reflecting their zone folding 
nature, the characteristics of Ai vary with n. 

The nature of the band gap in our superlattices was 
determined by the competition between the novel defect- 
derived flat state and the folded states of the Si(lll) 
layers in the conduction band. For small n, the band 


FIG. 3. (a)-(c) The PBE band structures are compared for 

the Si(lll)n/Si(SG) superlattices (n = 3-5) with the cubic- 
and hexagonal-stacking sequences of the Si(lll) layers. 


gap of the non-defective region becomes large due to the 
quantum confinement effect. Thus, the defect-derived 
fiat band at F is lower than all the folded states, ex¬ 
hibiting a direct band gap behavior up to n = 5. The 
planar-averaged charge densities clearly show that CBM 
is largely confined in the defective region, while VBM 
is mainly derived from the Si(lll) layers (Fig. |^. For 
large n, as the confinement effect is reduced, the folded 
Ai states move down below the defect-derived flat band, 
resulting in a direct-to-quasidirect band gap transition 
around n = 6 (Fig. [^. 

The novel defect-derived band plays an important role 
in strong optical transitions near the threshold energy. 
The squares of the dipole matrix elements, |(/|^^)P, for 
the direct transition at F were calculated to be 0.173, 
0.141, and 0.045 Uq ^ in atomic units for n = 3-5, re¬ 
spectively, where uq is the Bohr radius. These values are 
higher than 0.03 Uq ^ obtained from the dipole-allowed 
direct band gap of a specially designed Si/Ge super¬ 
structure [11], indicating that the optical transition was 
greatly enhanced in our case. Such strong dipole-allowed 
transitions are partly attributed to the large overlap of 
the band edge states at the interface layers (Figs. and 
1^. In addition, the hybridization of the conduction band 
states is another critical factor for the dipole-allowed 
transition. In c-Si, the band edge states at F are charac¬ 
terized by bonding and antibonding p orbitals and the op- 
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FIG. 4. The PBE band (left panel) and atomic (right panel) 
structures are compared for (a) the Si(lll)n=5/Si(SC) super¬ 
lattice with the cubic-stacking sequence of the Si(lll) layers, 
(b) a slab geometry consisting of hve Si(lll) layers and a vac¬ 
uum region, and (c) a slab geometry consisting of a defective 
layer sandwiched between two Si(lll) layers and a vacuum re¬ 
gion. Surface Si dangling bonds are passivated by hydrogen. 
In right panel, isosurfaces (yellow) of the charge densities of 
the lowest conduction bands at the F point are drawn. 


tical transition at the direct gap is dip ole-allowed. In our 
superlattices with direct band gaps, we examined the or¬ 
bital characteristics by projecting the wave function onto 
the atom-centered spherical harmonics within a sphere 
of 1.5 A radius. In flanking Si(lll) layers, tetrahedral 
bonds, especially around the interface Si atoms bonded 
to the SCs, are distorted from their ideal values, with 
the bond angles of 90.4° and 128.5° (Fig. [^. As a conse¬ 
quence, the p orbital character is significantly enhanced 
for the CBM state, as in the dip ole-allowed transition at 
the direct gap of c-Si. The calculated absorption coef- 



FIG. 5. The lowest conduction band energies at F, Ai, and 
A 2,3 are plotted as a function of n and compared with that 
of c-Si (horizontal dotted line). 



FIG. 6. For the Si(lll)n=5/Si(SG) superlattice, (a) isosur¬ 
faces (1.32 X 10^ electrons/A^) of the charge densities of VBM 
(blue) and GBM (red), (b) local density of states (in units of 
electrons/eV/A^) averaged over the xy plane, and (c) planar- 
averaged charge densities of VBM and GBM are plotted along 
the superlattice direction (z-axis). Black dotted lines denote 
the position of interface Si(lll) layers and white dotted lines 
represent the approximate positions of the band edge states 
in the middle of the defective and non-defective regions. 
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FIG. 7 . For various Si(lll)n/Si(SC) superlattices with the 
cubic- and hexagonal-stacking sequences of the Si(lll) layers, 
(a) the calculated absorption coefficients are compared with 
the experimentally measured values (from Refs. [30H33]) for 
GaAs, GdTe, GuInS 2 , amorphous Si (a-Si), and c-Si, and (b) 
the spectroscopic limited maximum efficiency (SLME) is plot¬ 
ted as a function of film thickness L. 


ficients of our superlattices are comparable to those of 
direct band gap semiconductors, such as GaAs, CdTe, 
and CuInS 2 [Fig. [^a)], which are known as good pho¬ 
tovoltaic materials. The spectroscopic limited maximum 
efficiency [34] for the sample thickness of L = 2.0 jam was 
estimated to be in the range of 27%-31% for n = 3-7 
[Fig.^b)], indicating that the optical absorption prop¬ 
erties were excellent even for the quasidirect band gap 
superlattices. 

Note that our superlattices are both energetically and 
dynamically stable at the ambient condition. Since the 
Si(lll) layers are structurally almost identical to c-Si, 
their excess energies are quite low, ranging 0.013-0.042 
eV per atom for n = 3-10 (Table [^. These energies 
are lower by an order of magnitude than those (0.1-0.3 
eV per atom) of the previously predicted Si allotropes 
with direct and quasidirect band gaps dMH]. For our 
direct band gap superlattices, we found no imaginary 
phonon modes in the phonon spectra (Fig. |^. In ad¬ 
dition, we examined their thermal stability by perform¬ 
ing first-principles molecular dynamics (MD) simulations 
and confirmed that they were thermally stable up to 100 




FIG. 8. For the Si(lll)n/Si(SG) superlattices with the cubic- 
stacking sequence of the Si(lll) layers, the calculated phonon 
spectra are shown forn = 1, 3, 4, and 5. The 2x2x3, 3x2x2, 
3x2x2, and 3x2x2 supercells are used for n = 1, 3, 4, 
and 5, respectively. 


ps at the high temperature of 1100 K (Fig. |^. 

Finally, we examined the possibility of creating an SC 
layer, which can lead to the eventual realization of our su¬ 
perlattice structures. We performed first-principles MD 
simulations for wafer bonding between two Si(lll) 2x1 
surfaces. For the initial configuration, where one surface 
had Pandey 7r-bonded chains [35] and the other SCs [16] , 
a 4x2 lateral supercell was chosen, with a vacuum region 
of 15.0 A inserted between the two surfaces passivated by 
hydrogen. It was found that a defective layer with the 
SCs was formed after about 3 ps at 1100 K [Fig. [T^a)]. 
Here, the 2x1 reconstruction of the SCs may not be pos¬ 
sible, since it is metastable with respect to the Pandey 
reconstruction model. On the other hand, divalent ad¬ 
sorbates, such as Ca, Sr, and Ba, can stabilize the SCs at 
half-monolayer coverage [36]. Thus, we considered one 
of the 2x1 surfaces with Ca atoms adsorbed at hollow 
surface sites [Fig. [Tol^b)] and found that the same defec¬ 
tive layer was formed after approximately 7 ps. In the 
final configuration, the Ca adsorbates resided along the 
open channels of eight-membered rings, as in the case 
of the Eu 4 Ga 8 Gei 6 -type structure found in CaSie, SrSie, 
BaSie, [37] and recently in NaSie [28l[38]. Once the de¬ 
fective layer is formed, the adsorbed Ca atoms can be 
removed via a diffusion process along the channel, as ex¬ 
perimentally observed in the case of NaSie [28] . Thus, the 
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FIG. 9. For the Si(lll)n/Si(SC) superlattices with the cubic- 
stacking sequence of the Si(lll) layers, the thermal stability 
was examined by performing first-principles MD simulations 
at 1100 K, with choosing the 2x2x3, 3x2x1,3x2x1, 
and 3x2x1 supercells for n = 1, 3, 4, and 5, respectively. 



FIG. 10. (a) For the initial configuration where one surface 

has Pandey 7r-bonded chains (blue circles) and the other has 
Seiwatz chains (orange circles), an intermediate conhguration 
after 2 ps and the final configuration after 3 ps at 1100 K 
are shown during molecular dynamics simulations for wafer 
bonding, (b) For the initial configuration where Ga atoms 
are adsorbed between Seiwatz chains on one surface, an inter¬ 
mediate configuration after 3 ps and the hnal conhguration 
after 7 ps at 1100 K are shown. In the hnal state, the Ga 
atoms are aligned along the open channels of eight-membered 
rings. 


wafer bonding between the clean and divalent-adsorbed 
Si(lll) surfaces can serve as a promising technique for 
the synthesis of superlattices containing the defective SC 
layer. 


IV. CONCLUSIONS 

Using a computational search method, we have pre¬ 
dicted low-energy pure-Si superlattice structures with 
dip ole-allowed direct band gaps, which can serve as 
promising materials for solar-cell applications. In the su¬ 
perlattice structure, the defective layers containing Sei¬ 
watz chains are intercalated in the diamond lattice of 
bulk Si. A direct-to-indirect band gap transition occurs 
as the bulk-like Si portion increases, which is understood 
in terms of a defect-originated conduction band and zone 
folding with quantum conhnement effects on the compet¬ 
ing conduction band of bulk-like Si layers. The method 
proposed for the synthesis of the defective layer through 
wafer bonding can be followed up by experiments and 
could herald a new optoelectronic era. 
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